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The behaviour of a miniature calorimetric sensor, which is under consideration for catheter- 
based coronary artery flow assessment, is investigated in both steady and pulsatile tube 
flow. The sensor is composed of a heating element operated at constant power, and two 
thermopiles that measure flow-induced temperature differences over the sensor surface. 

An analytical sensor model is developed, which includes axial heat conduction in 
the fluid and a simple representation of the solid wall, assuming a quasi-steady sen- 
sor response to the pulsatile flow. To reduce the mathematical problem, described by 
a two-dimensional advection-diffusion equation, a spectral method is applied. A Fourier 
transform is then used to solve the resulting set of ordinary differential equations and an 
analytical expression for the fluid temperature is found. To validate the analytical model, 
experiments with the sensor mounted in a tube have been performed in steady and pul- 
satile water flow with various amplitudes and Strouhal numbers. Experimental results 
are generally in good agreement with theory and show a quasi-steady sensor response in 
the coronary flow regime. The model can therefore be used to optimize the sensor design 
for coronary flow assessment. 



1. Introduction 

Flow sensors based on forced convective heat-transfe r, such as hot-film anemometers , 
can be used for th e assessment of art erial blood flow (Seed &: Woodlll97nl : IClarklll974l : 



Nerem et In a recent study, iTonino et all (|2009t ) showed that if the treatment 



of patients with coronary artery disease is based on an indirect measure for coronary flow 
(derived from coronary pressure measurements), the clinical outcome improves signifi- 
cantly. Clearly, direct flow assessment by miniature sensors that can be introduced into 
the coror iary arteries would prov ide even more information about the condition of these 
arteries ( van 't Veer et al. l2009l) . In this study, we aim to characterize the behaviour 



of such a miniature convective heat-transfer sensor in steady and pulsatile tube flow, 
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Figure 1. Schematic view of the calorimetric sensor mounted on the inside of the tube wall 
(not drawn to scale); heater in black, thermopiles in grey, r representing the radial, x the axial, 
and cj> the circumferential direction, (a) Side-view showing the placement of the sensor elements 
(heater and thermopiles measuring Th — Tf and Td — Tu) with respect to the fluid flow, (b) 
Cross-section showing the positioning of the flexible device inside the tube. 



through both an analytical and an experimental approach. The sensor is based on a 
calorimetric flow measurement principle: it consists of a small aluminium heating ele- 
ment of width 6=140 |^m, operated at constant power, and two polysilicon thermopiles 
that measure flow-induced temperature differences over the sensor surface. These sensor 
elements are embedded in a flexible polyimide substrate having a thickness of 10 |im; see 
figure [l](a). In order to use it for coronary flow assessment, the flexible device is bent 
around a catheter guide wire, which can be inserted into the coronary arteries. In our 
characterization study however, the device is mounted at the inner wall of a tube, to be 
able to subject it to a well-defined flow regime. The length I of the device is equal to 
approximately half the circumference of the inner tube wall; see figure [T]( 6). The tem- 
perature difference between two positions 100 |J.m downstream and 100 |J.m upstream 
from the heater centre (Td — Tu) is measured, as well as the heater temperature Th with 
respect to the ambient fluid temperature Tf far upstream (2000 ixm from the heater cen- 
tre). In absence of flow, heat transfer from the sensor to the fluid occurs solely through 
conduction, resulting in a symmetric temperature distribution over the sensor surface. 
If a certain fluid flow exists, the advective heat-transfer leads to an asymmetric tem- 
perature distributi on. The resulting temperature differences are a measure for the flow 
( Elwenspoek 19991 ). Flow reversal will lead to a sign change in Td — Tu, and hence can 
be detected, which is an advantage compared to the conventional hot-film anemometers 
( van Oudheusden fc Huiising 19891 ). 

Two important dimensionless parameters that appear in the study of thermal sensors 
in time-dependent flow are the Peclet number Pe, and the Strouhal number Sr. Here, Pe 
is a measure for the importance of advective compared to conductive heat-transfer, Sr 
for the importance of unsteady compared to advective temperature variations. Formal 
expressions for Sr and Pe are given further on (see (|2.8|) and (|2.13[) . respectively). 

Many experimental and analytical studies of hot-film anemometers have been reported 

in the literature. Experimen t s with these k inds of probes have bee n per formed by, among 

others. lSeed fc Woodl(fl970l).|Clarkl ( ll97'j).lAckerberg et al\ IHzi), andlvan Steenhoven fc van de BeuckenI 



(|1991| ) . iLiepmann fc Skinner! (|1954[ ) derived a theoretical relation between the amount 



of convectiv e he at loss from a hot surface and the local steady wall-shear rate. IPedlevI 
( I972I . 19761 ) and Menendez fc Ramaprian (1985) extended this work to include unsteady, 
pulsatile fiows. 

For miniature flow sensors like the one presented here, the theory developed for hot-film 
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anemometers is not applicable. First of all, because in a thermal sensor with small dimen- 
sions, generally operated at small Pe, heat conduction in the flow dir ection cannot be ne- 
glecte d , as is done in t he usual boundary-layer approxi mation (see e.g. lLiepmann fc Skinneii 
1954 [) : |Pedlevl(ll972l):lMenendez fc Ramaprianl(lf985l )). In numerical studies. iTardu fc PhamI 



20051 ) and iRebav etall ^200'£ showed that this axial conduction has a considerable 



influence on the response of small hot-film gauges. [Ackerberg et al. (1978) derived an 
analytical solution f or the heat transfer from a finite strip for Pe — >■ 0, in steady flow. 
Ma fc Gernerl ( 1993 ) examined the leading and trailing edge of a micro-sensor in steady 
flow separately to obta hi an analytical solution for th e entire sensor surface, in analogy 
to the method used bv lSpringer fc Pedlev ( 1973 ) and Springer (1974). 



Second, most analytical studies consider the problem of a uniform surface tempera- 
ture on the heated element, while our sensor is ope rated at constan t po wer, which is 



better described by a heat-flux boundary condition. iLiu et al\ (|1994[ ) and iRebav et al 



(|2007( ) used a constant heat-flux boundary condition on the surface of a heated element 
embedded in an adiabatic wall. However, the thermopiles of our sensor consist of con- 
ductive polysilicon, and therefore heat will not only be transferred from the heater to 
the fluid directly, but also via the surrounding material. In that case, we end up with 
a conjugate heat-transfer problem, where the heat source is known , but the interface 
temperature and hea t flux to the fluid a re unknown; iTardu fc PhamI (|2005r ) studied this 
problem numerically. IStein et al. ( 2002 ) deri ved an analytical solution downstream of a 
flush-mounted heat source in steady flow. IColc (2008.) considered conjugate heat-transfer 
from a steady-periodic heated film, also in steady flow. To t he authors' knowledge, on ly a 
few analytical models specific for calorimetric sensors exist. iLammerink et al. ( 19931 ) de- 
scribed an experimental study and a relatively simple analytical model of such a sensor, 
but in steady flow and with a different sen sor geometry. Calorimet ric sensors on highly 
conductive silicon wafers are described by Ivan OudheusdenI (|l99ll ). also in steady flow 
and with very small on-sensor temperature differences compared to the sensor overheat. 
Ex perimental studies with calo rimetric flow sensors in steady flow have been reported 
bv ILammerink et al\ ( 1993 ) and lNguven fc Kiehnscherj ( 1995 ). However, no data on un- 



steady flow experiments with this type of sensor are available. 

Our analysis of the calorimetric flow sensor focuses on the derivation of a new analytical 
model for the temperature distribution in a pulsatile fluid flow over a small heated element 
operated at constant power. Experiments with the sensor in steady and pulsatile water 
flow with Strouhal numbers and amplitudes in the expected physiological flow range 
are carried out to verify our theoretical predictions. In Jj^jthc mathematical formulation 
of our problem is given in terms of a two-dimensional advection-difFusion equation. We 
circumvent the coupling of the heat-transfer problems in the fluid and the substrate by 
approximating the heat flux from the sensor to the fluid, and use this approximation as 
a boundary condition for the fluid compartment. The axial conduction term is retained, 
and therefore our solution holds for all Pe-values. Since the heat flux at the boundary 
is approximated by a continuous fun ction, the leading an d trailing edge of the heater 
do not have to be treated separately (|Ma fc Gerneil 119931 ). resulting in one solution for 
the complete domain. When applied in coronary flow, our sensor will be operated at 
small Strouhal numbers, therefore a quasi-steady sensor response to the pulsatile flow is 
assumed. As described in ^ a spectral method is applied to reduce the mathematical 
problem to one dimension. Then, a Fourier transform is used to solve the resulting set 
of ordinary differential equations. The experimental technique is described in 2] In © 
the experimental results are compared to the theoretical predictions, and found to be in 
good agreement. The model developed not only leads to theoretical understanding of the 
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operating principle of the sensor, but it can also be used to optimize the sensor design, 
as is demonstrated in ^ 



2. Mathematical problem formulation 

In order to formulate an analytical model for the sensor in a pulsatile tube flow (see 
figure [1]), a cylindrical coordinate system (r, 0, a;) is adopted, where the main flow is in 
axial or x-dircction, r is the radial, and (j) the circumferential coordinate. The origin 
of this system is chosen such that x ~ Q &i the heater centre. The pulsatile fluid flow 
is assumed to be fully developed, and, because the temperature difference between the 
heater and the oncoming fluid is relatively small, temperature independent. The typical 
buoyancy-driven radial velocity can be estimated from the momentum equation in the 
radial direction using the Boussinesq approximation. For our configuration, the ratio of 
radial to axial velocity is of order 10^^, hence free convection can be neglected. 

The basic problem is thus reduced to that of finding the temperature distribution 
T(x,r, (f>,t), with t the time, in a prescribed pulsatile fluid flow in a tube of radius i?, 
which is heated by a time-constant prescribed heat influx in a small region of length I 
and width b around the tube wall; the remaining part of the wall is thermally insulated. 
The equation governing the temperature distribution in the fluid is the thermal energy 
equation in the tube {x £ M.,0 < r < R, — tt < (f> < tt) 

'd^T 1 d ( dT\ 1 d^Tl 



dT , ,dT 



Q^2 ' r dr \ dr 



(2.1) 



together with the boundary condition at the tube wall, r = R, describing the prescribed 
heat influx, 

fc_(a:,i?,^,i) =q(x), if ^ <0< — , (2.2) 
= 0, ifH>i^. 

Here, u is the velocity in z-direction, a the thermal diffusivity, and k the thermal con- 
ductivity. Heat influx q is in Wm^^, such that the power supplied to the heater in W is 
given by 



Q = l q{x)dx. (2.3) 



Considering the case that u > 0, i.e. the fluid is flowing in the positive a;-direction, we 
state that T must go to To, the initial fluid temperature, for a; — > — oo, but that T for 
X — > -|-oo must tend to a value T^^ > Tg for a quasi-steady solution to exist. In that case 
the total heat-transfer rate Q into the fluid is balanced by the advective heat outflow in 
positive a;-direction, equal to pcD{Too — To), with p the density, c the specific heat, and 
D = 2tt u{r, t)rdr, the total volumetric flow rate at time t. 

To emphasize the effect of the two different length scales that arise in the problem, i.e. 
heater width b and tube radius R, we introduce the following dimensionless variables 

X r u ^ t - T — Tn 

' = V "=y' ^ = ^' ^'-'^ 

with V the typical axial velocity, tc the time scale for temperature variations, and Tc the 
typical temperature scale. The characteristic parameter values can be found in table [TJ 
appropriate choices for tc and Tc are explained below. Note that the heater width is used 
as the characteristic length scale in x-direction, implying that we will look for changes 
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Parameter 


value 


unit 


description 


To 


20 


°C 


outer flow temperature 


To 


11.7 


°c 


temperature scale 


Q 


80 


mW 


heater power 


Xh 





(i.m 


heater centre 


Xd 


100 


fi.m 


position where Td is measured 


3ju 


-100 


(i.m 


position where Tu is measured 


Xf 


-2000 


|i.m 


position where T; is measured 


b 


140 


|i.m 


heater width 


I 


7000 


|i.m 


heater length 


cr 


70 


|i.m 


standard deviation of the assumed 
boundary heat-flux distribution 


R 


2.5 


mm 


inner tube radius 


V 


0.1 


ms~^ 


typical axial velocity 


So 


115 


s-i 


mean wall-shear rate 


a 


1.44 ■ 10"^ 


m^s-i 


thermal diffusivity^ 


k 


0.606 


Wm~^K~ 


^ thermal conductivity^ 


V 


1 • 10"" 


m^s-i 


kinematic viscosity :j: 


w 


27r 


rad s~^ 


angular frequency 




X (see 


llncroDera 


et aZ.ll2007l. d. 860) 



Table 1. The parameter values used in the analytical model, based on the experimental 

set-up. 



in temperature T in the direct axial vicinity of the heater, which is where and are 
measured. By substituting (j2.4|) into we obtain 



h df . df _ a d^f 

_._+u{r,t)—-— — 



1 d^f 




f2 Q^2 



(2.5) 



with e = \/ {ab/V)/R = 0.006 ^ 1. Since, in (|2.5p . the small number appears in front 
of the highest derivative with respect to f, one can expect a boundary layer to develop 
at the tube wall; i.e. at f = 1. The outer solution at leading order, with e = 0, is the 
trivial solution T = 0. The temperature problem is thus confined to a small region close 
to the sensor surface: the thermal boundary layer of thickness 5t- 

Our sensor measures the temperature difference Td — T„ at distances in the order of 
magnitude of b up- and downstream of the heater centre for |(/)| < 1/2R. Hence, in the 
region of interest for our sensor |x| = 0{b) and 5t <C R- Within this region, the problem 
is independent of the (/i-coordinate. The characteristic length-scale for conduction in the 
0-direction, heater length / = 0{R), is much larger than the length scale for conduction 
in the axial direction, heater width 6. For our sensor b/l = 0.02; hence conduction in the 
<?!>-direction and edge effects occurring at (j) ~ ±l/2R can be neglected. 

Since the thermal boundary layer thickness 6t is much smaller than the tube radius, 
the tube wall in a 6-environment of the heater can be considered flat. We therefore 
adopt a spatial rectilinear coordinate system (x, y), where x is the surface coordinate in 
the flow direction, and y is the stretched coordinate normal to the surface, defined as 
y — (i?/(5T)(l — r); see figure [2] As a further approximation, we confine the domain for 
the inner solution to a strip of finite height h. At the upper boundary of the strip, y — h 
(with y = dry), we then require that T = 0, to match the inner solution to the outer one; 
see figure [2l How to choose h such that the solution in a &-environment of the heater. 
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y=h: upper boundary 





u(y,t)=S(t)y / y 






./ 


T=T„for x-> ±oo 






y=0: sensor surface 


V Y T 

- 1 iX^ 1 


I A A A _ * 


^ - -1 





Figure 2. Scheme of the problem geometry. 



where the thermal boundary layer is still thin, is not influenced by the finite size of the 
domain in y-direction will be explained further on in this section; see (|2.1ip . 

Further downstream (for x > b, hence outside our region of interest) the thermal 
boundary layer widens, due to radial conduction. Both curvature and (/)-dependence will 
enter the problem again, while axial conduction will become negligible. Even further 
downstream {x > R/b), the fluid temperature will become uniform in each cross-section, 
with T — > Too ^ Tq + Q/ pcD. It is therefore important to note that, given the simplifica- 
tions described above, our method will only yield the correct solution in an 6-environmcnt 
of the heater, i.e. the region of interest for our sensor. 

As a further approximation, we assume that the wall-shear rate is the only flow param- 
eter that influences the heat transfer from the sensor surface, implying that the velocity 



profil e may be approximated linearly throughout the thermal boundary layer (see lPedlev 



19721) . Since our domain is now restricted to a strip of finite height h, the linearization of 



the velocity profile is valid throughout the complete domain. This approximation requires 
the Stokes layer thickness Ss to be much larger than the thermal boundary layer thickness 
^7^. In that case the velocity u within the thermal boundary layer can be approximated 

by (y = - 7- = Stv) 



y = S{t)y, 

with S the wall-shear rate, which is, in a pulsatile tube flow, given by 



S{t) = Sq[1 + 13 sin {ujt)], 

with 5*0 the mean, cj — 27r/ the angular frequency, and /? the amplitude of the shear-rate 
oscillations. This implies that when /? > 1, backflow is involved. Since for coronary flow, 
the order of magnitude of (3 will be about one, the dimensionless shear rate S{t) = S{t)/ Sq 
is an 0(l)-function of t. Hence, we have 

u{y,t) = S{t)y = SoSrySii), 

yielding 

V = StSq, and u{y, i) = S{i)y. 

To find an expression for the thermal boundary layer thickness St, we write (|2.5|) in 
terms of y, u{y,t), and S{t) as 

b df ^,^,^df a d^f ab d^f 

^WVl^ = uo , ^^2 + .3c a..2 - (2-6) 



tcSoSx di ' ' dx bSodr dx^ S^Sq dy^ 



Characterization of a miniature calorimetric flow sensor 



7 



When advection in x- and conduction in y-direction are the two dominant effects, ab/ SoSj- 
must be of 0(1), and he nce the thermal boundary layer thickness is given by (see 



Liepmann fc Skinneij|l954l ) 



St = {ah/Sof^ . (2.7) 

The linearization of the velocity profile within the thermal boundary layer is allowed if 
5t 5s- The Stokes layer th ickness in a fully developed pulsatile tube flow is given by 
(|SchUchting fc Gerstenll2000[ p. 367) 



5s = {v/u^f^ , 

with V the kinematic viscosity. For a 1 Hz pulsatile water flow we get (5s = 4 • 10^^ m. 
The requirement 5t ^ 5s leads us to an estimate for the admissible shear rate: 

5o > ah/5% = 0.3 s'^ 

This requirement is amply satisfied, since our experiments are performed at a mean wall- 
shear rate of about 115 s~^. Furthermore, 5t and b are of the same order of magnitude 
in this range of shear rates, allowing the use of b as length scale in both x- and y- 
direction. This already indicates that the axial conduction term in (|2.6|) can certainly 
not be neglected within the region of interest for our sensor. This is further confirmed 
by the magnitude of the coefficient of x-conduction; a/bSoSx = 0.16 (see table[l]for the 
parameter values used). 

The magnitude of the dimensionless group in front of the unsteady term in (|2.6p . the 
Strouhal number Sr = b/tcSodx ~ ^/tcSo, if b/Sx ~ 1, depends on the choice of the 
characteristic time scale tc- The goal of this study is not to analyse start-up processes 
that occur when switching on the heater, but to describe the periodic variations in the 
sensor response. Therefore, the oscillation time, l/uj, is used as characteristic time scale; 
hence 

Sr = w/Sq. (2.8) 

Eventually, this sensor will be used for coronary flow measurements, where the estimated 
mean shear rate the sensor experiences when positioned on a guide wire is in the order 
of magnitude of 1000 s~^. Hence, Sr is generally small for our ultimate application 
(typically Sr < 0.1, assuming measurem ent of the fir st ten harmonics is sufficient for 
reconstruction of the coronary flow signal (|Milnoilll989l p. 157)). We therefore assume the 



fluid temperature distribution to be quasi-steady, thereby neglecting the unsteady term 
in the thermal energy equation. In the quasi-steady approximation time t represents a 
parameter rather than a variable. From here on, we therefore omit the explicit dependence 
on i (T = T{x,y)); in fact, the role of t is now taken over by the shear rate S {S = 
l + /3sin(t)). 

Another simplification is that heat loss through the insulating back of the tube in 
which the sensor is mounted, is neglected: all heat produced by the heater is assumed 
to be transferred to the fluid. Capacitive effects, which ma y cause the heat transfer to 
the fluid to vary in time, are also neglected. According to iTardu fc Phain (2005) this 



is reasonable, since the thermal diffusivitics of the sensor components are two orders of 
magnitude higher than that of water. Since the exact shape of the heat-flux distribution 
from the sensor substrate to the fluid depends on the temperature distribution in the fluid, 
this leads to a conjugate heat-transfer problem, which is hard to solve. We circumvent 
this coupling of the fluid and substrate temperatures by a much simpler approach: we 
approximate the shape of the heat-flux distribution from substrate to fluid and use this 
as a boundary condition for the fluid problem. If all the heat would be transferred from 
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the heater to the fluid directly, i.e. when there is perfect insulation outside the heater 
compartment, a rectangular-shaped heat-flux boundary condition would be most realistic. 
For our sensor, however, conduction of heat from the heater towards the other sensor 
components will smooth the rectangular-shape, leading to a more Gaussian-shaped heat- 
flux profile, with some deviations due to the asymmetric temperature distribution in the 
fluid. As a simple approximation of the real heat-flux boundary condition, we therefore 
use a Gaussian distribution with a standard deviation a equal to half the heater width, 
hence a = 6/2: 



q(x) = -k- 

dy 



with Xh the position of the heater centre, and Q/l the total amount of heat transferred 
from the sensor to the fluid, per unit of length in z-direction (in Wm^^), as given by 
((231) . 

The resulting dimensional thermal energy equation and boundary conditions describing 
the quasi-steady problem for the fluid temperature T = T{x, y) within the thermal 
boundary layer or strip are given by: 

dT fd^T d^T 

I = (T ^ -I ^ 



ox \ ox'^ oy'^ J 

r(±oo,y) = To, — (x,0) = - ^ e--^, T(x,h)^TQ. (2.10) 
ay klaV^n 

Figure[2]shows a schematic view of the resulting problem to be solved. Note that, although 
we are only interested in the solution close to the heater (for |a;| = 0{b)), we have, for 
mathematical ease, extended the domain in x-direction to infinity. The outflow boundary 
condition used, implies that in our model, all heat will eventually escape through the 
upper boundary y = h. Since, for the correct choice of h, this happens sufficiently far 
away from the heater, it does not influence the solution in a ^-environment of the heater. 
To ensure this, the upper boundary of the domain has to be located sufficiently far 
outside the thermal boundary layer for x = 0{b). Therefore, h is taken equal to n times 
the estimated thermal boundary layer thickness, where n = 1,2,3,.... Further on, it 
will be shown that for n > 4 the solution becomes independent of n. Since the thermal 
boundary layer thickness St according to (|2.7p . but with 5*0 replaced by S, depends on 
the actual wall-shear rate, h depends on S as well. This motivates us to choose h as (note 
that in the following we use o-\/2 instead of b as characteristic unit of length) 

Hence, in our quasi-steady approximation we solve the problem for each value of S 
separately, choosing the upper boundary accordingly. The advantage of this S'-dependcnt 
position of the upper boundary will become clear in the next section. We note that h 
can become large, i.e. larger than 6, for small values of \S\, specifically for \S\ < 0.3 s~^. 
Then, the thermal energy equation is no longer advection-, but diffusion-dominated and 
in that case St must be taken equal to b. However, for our problem this happens in a 
very short period of time (less than 0.1% of one period of S{t)) and it is therefore not 
relevant for our solution. 

We introduce a new scaling by using the dimensionless variables and parameters: 

x-xh . y T{x,y)-To 

X = y = — 1=1 J {x,y) = — , 

aV2 aV2 
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„ Q r h a , ^ 

h=^, ^ = 7rT^^ (2-12) 



fc/v^' aV2' 2ct2S" 

where the temperature scale Tc is based on the heat source term (|2.9p . We define the 
Peclet number as 

Pe.^, (2.13) 
a 

hence a=l/Pe. 

Omitting the tildes, the newly scaled system for T = T(x, y) reads 
dT /d^T d^T\ 

dT 2 

T(±oo,y) = 0, — (x,0) = -e-% T(x, /i) = 0, (2.14) 

where 

h = (2.15) 

with n still to be chosen. Hence, T{x,y) depends on only two parameters, a and n: 
T{x,y) ~ T{x,y]a,n). However, if n is taken sufficiently large, i.e. n > 4, then solution 
T in a 6-environment of the heater becomes independent of n, and a is the only parameter 
remaining. 



3. Analytical solution method 

To solve the system (|2.14p a spectral method is used, which reduces the partial differ- 
ential equation in (|2.14l) to a set of ordinary differential equations. To apply this method, 
we first make the boundary conditions homogeneous, by writing 

T{x,y)^{h~y)e-^" +T,{x,y), (3.1) 

leaving for Ti the equation 

dTi fd^Ti d'^T{ 

a \ 

dx \ dx'^ dy"^ 

with homogeneous boundary conditions and with 

R{x, y)^{h- y) e--" [2xy + a {4x^ - 2)] . (3.3) 
For the spectral method, we introduce the trial functions Vk{y), given by 

^ = -A^,, ^(0)=0, vu{h)=0. (3.4) 
dy Ay 

yielding 

(2k — l)7r 

wfc(y) = cos(Afcy), Afc = ^ , fc = l,2,.... (3.5) 

Here we see the advantage of truncating the infinite half-space to a strip of finite height. 
Next, we decompose T\ into a linear combination of the trial functions v-^ according to 

oo K 

Ti{x,y) = ^Ck{x)vk{y) ~ ^Ck{x)vk{y), (3.6) 

k=l k=l 

where in the last step we have truncated the series after K terms (as demonstrated in fJU 
iiT = 5 is more than sufficient for obtaining precise numerical results when the height of 
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the strip is chosen according to the actual wah-shear rate). Substituting p.ip and p. 6 
into ()2.14|) . we obtain 



K r 



E 

1=1 



y- a 

ax 



(3.7) 



Taking the inner product of (|3.7p with functions Vk{y), with the inner product of a 
function u with v defined as 

(u,u) EE / u{y)v{y)Ay, 
Jo 

we arrive at an equation for the array C, consisting of K elements Cfc, 



'^^d^-"2 d^ 



-AC\=R 



(3.8) 



with W a. K X JT- matrix with elements Wki, given by p.9a[) . A a. K x K diagonal matrix 
with elements A^^ ~ \\ and R a iiT-array with elements Rk, given by p.96p : 



Rkix) 



yvk{y)viiy)dy ^ / yvk{hy)vi{hy)dy, 



Rix,y)vk{y)dy. 



We introduce the Fourier transform of C{x) by 

1 f°° 

c(C) = ^ / C(a;)e-<-dx = ^{C;C}. 



(3.9a) 
(3.96) 

(3.10) 



By taking the Fourier transform of (|3.8p . after dividing it by ah/2, we obtain the algebraic 
equation for c: 



+ ^(^^W + A] c(C) - M(C)c(C) - r(C), 



with / the unity K x A'-matrix, and 



2h. 



M{C) = CI + iC—W + A, 

a 

.(C) = ^{|^i^;C}. 



(3.11) 

(3.12a) 
(3.12&) 



We can. using Mathematica 6 (Wolfram Research, Champaign, USA), invert the K x K- 
matrix M analytically, by which we find 



c(C) = M-i(C)r(C), 



(3.13) 



and by taking the inverse Fourier transform of this result, we obtain the solution for the 
array C{x) as 



1 /'°° 1 /'°° 

C(x) = -= / c(C)e<-dC=^ / M-\C)r{OS^dC 



(3.14) 



The latter integral is evaluated numerically using Mathematica 6. The temperature T is 
now determined by p.ip and p.6p with K = 5 and n = 4. 
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4. Experimental methods 

In the experimental set-up, the device was mounted to the inner wall of a tube with an 
inner diameter of 5 mm (sec figure[T](6)). The tube was made of PMMA, which is an in- 
sulating material, to prevent heat loss through the back of the device. The device covered 
about half of the tube perimeter. The polyimide foil including the sensor components 
has a thickness of only 10 i^m (see figure[T](a)), and since the very small step it causes in 
the tube wall is located about 2 mm away from the actual sensor components (on both 
sides) this does not significantly disturb the flow pattern near the sensor. In all experi- 
ments, the heater was supplied with a power of 80 mW by a voltage source (EST 150, 
Delta Elcktronika, Zicrikzee, The Netherlands). Two multimeters (DMM 2000, Kcithley 
Instruments Inc, Cleveland, USA) were used to register the output of the thermopiles 
that measure — T.^ and — Tf. 

To ensure fully developed flow over the sensor, the measurement section was located 
112 tube diameters from the tube entrance, which is, even at the highest Reynolds num- 
ber reached (« 500), well beyond the laminar entrance region. As test fluid, tap water 
at room temperature was used. Steady flow through the set-up was generated by a sta- 
tionary pump (Libel-Project, Alkmaar, The Netherlands). The amount of flow could be 
adjusted using a clamp. The oscillatory component was added to the mean flow by a pis- 
ton pump, driven by a computer-controlled motor (ETB32, Parker Hannifin, Offenburg, 
Germany). Downstream of the sensor, the flow was registered by an ultrasonic flow probe 
(4PSB transit time perivascular probe. Transonic Systems Inc, Ithaca, USA), which was 
used as a reference. The signals from the flow probe and the multimeters were recorded 
simultaneously and transferred to a computer via an acquisition board with a sampling 
frequency of 20 Hz. 

The output voltage of a thermopile is proportional to the temperature difference be- 
tween its ends via the Seebeck coefficient, which d epends on the composition of the 
thermocouple leads (see van Herwaarden et all [l989l ) . By scaling the stationary sensor 
response Th — Tf at S = 115 s~^ to the analytical value at this shear rate, we found the 
Seebeck coefficient for the thermocouples in our sensor to be 305 |xVK~^. This Seebeck 
coefficient, which is only a scaling value for the experimental data and docs not influence 
the shape of the responses, is used for all experimental results shown in fJS] 

Experiments under both steady and unsteady flow conditions were performed. In 
steady flow, the Peclet number (see (|2.13p ) was varied from (at zero flow; 5* = 
s~^) to 34 (at a flow of 368 ml min~^; S = 500 s~^). The wall-shear rate at the sen- 
sor surface was calculated from the flow measured by the ultrasonic probe assuming a 
Poiscuille velocity proflle. 

For the unsteady case the Strouhal number (see (|2.8p ) was varied from 0.01 to 0.1 by 
varying the oscillation frequency from 0.2 to 2 Hz, and the amplitude (/3) was varied 
between 0.8 and 1.2 (corresponding to the expected coronary flow regime), keeping the 
mean shear rate constant at about 115 s~^. In unsteady flow, the wall-she ar rate was 
deriv ed from the flow measurements assuming a Womersley velocity profile ([Womerslev 
19551) . 



5. Results and discussion 

The cosine-series used in (j3.6p converged quite rapidly: only five terms sufficed for 
an accurate approximation of the solution; see figure [3]( a). The rapid convergence is a 
consequence of the dependence of h on the actual wall-shear rate S; if h would have been 
fixed for all 5, it could become much larger than the boundary-layer thickness, leading 
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Figure 3. The influence of (a) the number of cosine-terms K on the solution Td — Tu and (b) 
parameter n on the solution T at y = 0, 5* = 100 s~^. 

to slow convergence of the cosine series. The parameter n (see (|2.15p ) was chosen such 
that the position of the boundary condition did not influence the solution at the sensor 
surface in a 6-environment of the heater; n = 4 was found to be large enough to ensure 
this, as demonstrated in figure [3](6). 

In figure |4](a), the theoretical temperature profiles over the sensor surface are de- 
picted for different shear rates (hence different Pe-values). At low shear rates (low Pe), 
the temperature distribution is more symmetric with respect to the heater centre, since 
in that case conduction is dominating the heat-transfer process. As the shear rate in- 
creases the temperature distribution becomes asymmetric, because more heat is advccted 
downstream, while the overall sensor temperature decreases because of the augmented 
advectivc cooling. 

The experimental data obtained in steady flow arc plotted together with the theoreti- 
cally predicted sensor output in figure |3](6-(i). The results of two separate experiments 
are shown to give an indication of the data spreading, where each data point represents 
the average result of 20 s of measurement with a sampling frequency of 20 Hz (i.e. 400 
samples). Both the experimental and analytical curves show a steep decline in the rel- 
ative heater temperature Th — Tf at low shear rates, and a more gradual one at higher 
shear rates. Sensor output — Tu is in both model and experiment characterized by a 
steep increase at low shear rates, f ollowed by a maximum a nd a decline; see figure [41(c) 



These features were also found by Lammerink et al. (1993) and Nguven fc Kiehnscherj 



([1995!). With increasing shear rate, the {Td — rii)-tcmperature difference rises because 
of augmented advection of heat in downstream direction. At the same time the overall 
sensor temperature decreases (see also figure |3]( a)), hence, a maximum in Td — Tu is 
observed. Apart from describing very well the qualitative sensor response, the analytical 
model also predicts the quantitative data with adequate precision. Maximum deviations 
between model and experiment ranged from 5% for the (T/j— T/)-signal to 27% for T^— T„. 
Although also measurement inaccuracies may play a role, the discrepancy between the- 
ory and experiment is most likely due to the simplified modelling of the substrate: the 
Gaussian heat-flux distribution is only a rough approximation, since the heat transfer 
from the substrate to the fluid will be larger upstream than downstream, due to the hot 
thermal wake. Furthermore, the influence of conduction in the substrate decreases with 



increasing wall shear-rate (see lTardu fc Phamll2005l ) 



To obtain an invertible relation between the sensor output and the wall-shear rate or 
Peclet number, the ratio of thermopile outputs, (Td — Tu)/{Th — Tf), can be used; see 
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shear rate (s^ ) 



(b) 25 -| 
20- 




100 200 300 400 
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Figure 4. Analytical ( — ) and experimental (• • •, ★★★) results in steady flow, (a) Analytical 
temperature profiles at the sensor surface at S =10, 60, 110, 160, 210 s~^, (b) response of the 
thermopile measuring Th — Tf, (c) Td — Tu, and (d) the ratio of thermopile outputs. 



figure 2] (rf). Note that this curve is independent of the thermopile calibration, because the 
Seebeck coefBcient is equal for both thermopiles and vanishes when the ratio of outputs is 
used. From figure 0] it appears that the sensor is most sensitive to lower shear rates. The 
performance of the sensor at higher shear rates, important for the eventual application of 
the sensor in coronary fiow, can be improved by decreasing the heater width 6, thereby 
reducing Pe. When the heater width = 2(t is decreased, the distance to the heater 
centre of the thermopile measuring Td — must be reduced by an equal amount, to 
keep the same relative positions. The theoretical results for decreasing the heater width 
by 50% and 25% are shown in figure [S] A smaller heater leads to a shift in the maximum 
temperature difference — T„ , resulting in a more linear relation between the shear rate 
and the ratio of thermocouple outputs, with a lower sensitivity for lower, and a higher 
sensitivity for higher shear rates compared to the original response. From figure a) we 
conclude that the effective heater width has a large influence on the (T^j — ru)-signal. 
This could also be an explanation for the discrepancy between theory and experiment 
shown in figure |4](c); if the effective heater width in the experiment is somewhat smaller 
than the theoretically used value, this will shift the maximum in the {Td — rti)-curve to 
higher shear rates. The difficulty here is that the effective heater width will depend on the 
actual wall-shear rate (i.e. the relative influence of conduction in the substrate), making 



14 



H. Gelderblom et al. 




Figure 5. Theoretical results with the original heater width ( — ), 50% (• ■ •), and 25% (- • -) of 
the original width for (a) Ta - T„ and (b) (T^ - Tu)/(Th - Tj). 



b a function of S. The actual effective heater width can therefore only be obtained by 
solving the conjugate heat-transfer problem. 

The sensor response to unsteady flow was investigated experimentally by varying the 
oscillation frequency, and thereby the Strouhal number, and amplitude in the estimated 
physiological regime. At each amplitude and frequency, the dynamic sensor response was 
measured during at least 5 flow cycles; here two periods of each signal are shown. In figure 
[H] the experimental results for non-reversing shear rates at four Sr- values are depicted, 
together with the quasi-steady analytical solution. The shear-rate signals calculated from 
the flow measurements are aligned, to ensure that the phase differences observed in the 
thermopile signals are due to thermal unsteady effects, and not to phase differences 
between the flow and the wall-shear rate. Owing to limitations of the pump, the 2 Hz 
{Sr~0.1) flow signal was not purely sinusoidal, which also shows in the sensor response. 

We observe a phase shift and decrease in amplitude with increasing Strouhal number in 
the (T^ — T/)-signal. The (T^ — Tu)-thermopile output also shows this phase shift, together 
with a slight change in the signal shape. As Sr increases, the wall-shear rate oscillations 
become too fast for the thermal boundary layer to react instantaneously, and the sensor 
response starts to deviate from its quasi-steady behaviour, and hence from the analytical 
solution. The deviation between the signals with Sr — 0.01 and Sr = 0.1 is larger (17% 
in Th — Tf) during minimum shear rate, when unsteady effects are most important, than 
during maximum shear rate (6%), when advection dominates. Not only during minimum 
shear rate, but during the complete deceleration phase the spreading between the different 
iSr-curves is somewhat larger. To investigate whether flow instabilities in the deceleration 
phase can explain these deviations, a frequency analysis of the signals was performed. 
No coherent structures with a fixed frequency were found, so the origin of the deviations 
remains unclear. Nevertheless, the quasi-steady analytical solution appears to describe 
the sensor response quite well in the complete experimental range, up to Sr = 0.1, 
with again larger quantitative differences in the (T^ — T„)-signal than the {Th — Tf)- 
signal. In the coronary flow regime, with Strouhal numbers of about 0.01 for the first 
harmonic, a quasi-steady sensor response is therefore expected. In thei r studies with hot- 
film anemometers and electrochem i cal wa ll-shear probes, respectively, IClarkI (119741 ) and 
van Steenhoven &: van de BeuckenI ( 1991 ) found the quasi-steady regime to hold for Sr 



up to 0.2. 




Figure 6. Results in unsteady flow for jS — 0.8, with r being the period of a flow cycle; analytical 

( ) and experimental (•••) curves for Sr— 0.01 (blue), Sr^O.OS (red), S'r=0.06 (green), Sr=0.1 

(magenta), (a) Shear rate at the sensor surface obtained from the Womersley approximation of 
the measured flow, (b) response of thermopile measuring Th — Tf, (c) Td — Tu, and (d) the ratio 
of thermopile outputs. 

For /3 — 1.2, larger deviations between the sensor response for Sr = 0.01 and Sr = 0.1 
have been found during the reversal period and the deceleration phase; see figure [71 
A sign change in Td — T^, indicating shear-rate reversal, was clearly observed for the 
two lowest 5'r-values. During the reversal period, hot fluid from the thermal wake is 
carried back over the sensor, which is not taken into account in the analytical model and 
leads to further deviations from the quasi-steady response. As the shear rate approaches 
zero, the heat is carried away from the heater only very slowly, leading to large heater 
temperatures in the quasi-steady analytical solution, while Td — Tu, and therefore also 
{Td — Tu)/(Th — Tf), tend to zero, due to the symmetric influence of conduction. In the 
experimental data such large relative heater temperatures are never reached because it 
takes time to heat the fluid, due to its finite thermal diffusivity. Hence, only during the 
very short period of time where S is close to zero, larger deviations from the quasi-steady 
solution are observed in the {Th — r/)-signals. 



6. Conclusions 

An analytical model describing the response of a miniature calorimetric sensor to both 
steady and pulsatile tube fiow is developed. In experiments the sensor is subjected to a 
flow in the expected physiological range to verify the theoretical predictions. Steady flow 
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Figure 7. See figure[6l with /3 = 1.2. 



analytical and experimental results are in good agreement for the complete range of Peclet 
numbers studied. Hence, our two-dimensional model with the wall-shear rate at the sensor 
surface as the only flow parameter is sufficient for examining the steady sensor behaviour. 
Only a simplified model of the substrate in which the sensor is embedded was taken into 
account, by means of a heat-flux boundary condition. A conjugate approach will lead to a 
more accurate quantitative prediction of the temperature differences measured, however, 
our model has the advantage of a simple representation of the substrate, and still leads 
to an acceptable description of the sensor response. 

The quasi-steady analytical model predicts the sensor behaviour in non-reversing pul- 
satile flow with Strouhal numbers up to 0.1 quite well. Based on the experimental results, 
we conclude that the sensor response to coronary flow will be quasi-steady, except during 
the (short) periods of shear-rate reversal. The analytical model can therefore be used to 
optimize the sensor design for coronary flow measurements, as demonstrated in figure [5l 

This research was financially supported by the Dutch Technology Foundation STW, 
project SmartSiP 10046, Philips Research, and St. Jude Medical. 
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